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Abstract 

The diffusion limit of the linear Boltzmann equation with a strong magnetic 
£S) ■ field is performed. The giration period of particles around the magnetic field is 

assumed to be much smaller than the collision relaxation time which is supposed to 
be much smaller than the macroscopic time. The limiting equation is shown to be 
a diffusion equation in the parallel direction while in the orthogonal direction, the 
guiding center motion is obtained. The diffusion constant in the parallel direction 
is obtained through the study of a new collision operator obtained by averaging the 
original one. Moreover, a correction to the guiding center motion is derived. 
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1 Introduction 



The motion of charged particles under the action of a strong magnetic field is an impor- 
Q\ tant phenomena encountered in artificial and natural plasmas (Tokamakas, ionospheric 

plasmas, etc). In the presence of constant magnetic field B, a charged particle (electron 
in our case) of mass m and charge q with velocity v follow a helical trajectory around 
magnetic field line. The gyration radius of the particles, called Larmor radius Tl, is in- 

^^ ■ TflV 

versely proportional to the magnitude of the magnetic field (rz = )• This means that 

qB 

when the magnetic field becomes large, the particles get trapped along the direction of 
B. In addition, when an electric field E is applied, the particles experience a drift of the 
instantaneous centers of Larmor circles, called in the literature the guiding center motion, 
in the direction perpendicular to both the electric and magnetic fields. The velocity of 
this drift is given by the following formula 

Ex B 
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Direct simulations of Vlasov or Boltzmann equations in the presence of such large magnetic 
fields requires the numerical resolutions of the small position and time scales induced by 
the gyration along the magnetic field. Hence, the question of deriving approximate models, 
numerically less expensive, is of great importance. Various models have been developed 
to this aim like the "guiding-center" and the gyrokinetic models. The guiding center 
approximation consists in averaging the motion over the gyroperiod when supposing B 
goes to +00 (t*£ — > 0) [26] and the gyrokinetic approximation is a generalization of the 
guiding center one in the case of kinetic equations. For a complete physical review, we refer 
for instance to [27J 1151 EH1 EI] • Mathematically speaking, homogenization techniques were 
used in [191 ED] by E. Frenod and E. Sonnendriicker to justify the high magnetic field limit 
of Vlasov and Vlasov- Poisson systems. It is proven in [19] that the tridimensional guiding 
center approximation leads to a one-dimensional kinetic model in the direction of the 
magnetic field. The drift phenomenon perpendicular to the magnetic field is rigourously 
obtained in [20] for the two dimensional Vlasov equation on a sufficiently long time scale. 
Other asymptotic regimes of the Vlasov equation with strong magnetic fields are studied 
[2T1 [TBI E21 HEl [T7J. The gyrokinetic approximation of the Vlasov-Poisson system has also 
been intensively studied by F. Golse and L. Saint-Raymond for the two and tridimensional 
motions in various asymptotic regimes [2H [331 E51 EI]- 

In all the above references, the transport is assumed to be ballistic, which means 
that particles do not suffer any collision during their motion. We are here interested in 
regimes were collisions are important and we propose to study the interplay between the 
fluid limit induced by collisions and the high frequency gyrations around the magnetic 
field. To simplify the study, we shall only consider collisions with a thermal bath at a given 
temperature. In this case, the diffusion limit of the obtained Boltzmann equation leads 
to the so-called drift diffusion equation. The question of deriving diffusions equations 
from the Boltzmann equation has been widely studied during the last three decades in 
the context of radiative transfer [51 El [351 ED E] , in semiconductors and plasmas [321 [231 USB 
El [31 HI [HI [H [TO] and many other contexts that we do not mention here. In the presence 
of a magnetic field, an important quantity to be considered is the gyroperiod measuring 
the time that a particle of mass m with charge q and submitted to the constant magnetic 
field B makes a 2ir rotation around this field (T c = 2imi/qB). When the gyroperiod 
is much larger than the relaxation time, then the magnetic field effect disappear in the 
diffusion limit as will be seen later on. When T c is of the same order of magnitude as 
the relaxation time, the diffusion matrix has an antisymmetric component, generated 
by the magnetic field (this has been proven for collision with walls by P. Degond & al 
[HI [UJ [TJ] ) . Formal results for binary gas mixtures can also be found in [291 H21 H31 [9] . We 
shall consider in this paper the situation where the gyroperiod is much smaller than the 
relaxation time. In that case, one expects that the diffusion matrix is obtained through 
the analysis of an averaged collision operator. Also, the high frequency oscillations only 
occur in the direction perpendicular to the magnetic field while parallel velocity stay 
unaffected. Therefore, the behaviours of the solutions in these directions are expected to 
be different. We shall show in this paper that the transport is diffusive along the magnetic 
field, while it is dominated by the guiding center drift in the orthogonal direction. The 
diffusion coefficient in the parallel direction will be obtained by averaging the collision 
cross section around the magnetic field. 

The outline of the paper is as follows. In the next section, we introduce the seal- 



ing, notations and state the main results. The proofs rely on the study of the operator 
accounting for collisions and gyration around the magnetic field. The analysis of this op- 
erator is done in Section 3. Section 4 is devoted to the asymptotics of this operator when 
the scaled gyration period tends to zero. It is shown that this limit is well described by a 
collision operator with a cross section averaged around the magnetic field. This result is 
then used in Section 5 in order to prove the main results of the paper. Some concluding 
remarks are listed in the last section of the paper. 

2 Setting of the problem and main results 

The problem we are interested in is a singular perturbation of the three dimensional 
Boltzmann equation with a constant large magnetic field. The position variable is denoted 
by r = (x,y,z), the velocity variable is denoted v = (v x ,v y ,v z ). The magnetic field is 
assumed to be constant and parallel to the z axis. We shall denote by rj_ = (x, y) and 
v ± — {v x ,v y ), the orthogonal variables. The distribution function is a solution of the 
following singularly perturbed Boltzmann equation 

df^ + %f^ + TlT + (v x e,) = OCT) 

dt E El] E 2 Tff V ' E 2 (2.1) 

/*"(* = 0) = / (r,v). 

where 71 and T z are the respectively perpendicular and parallel parts of the transport 
operator 

71 = vi ■ V rx + E ± • V v± , T z = v z 3 z + E z d Vz . 

The collision operator is given by 

W)(v) = / cr(v,vO[M(v)/(vO -M(V)f(v)]dv'. (2.2) 

Jm.3 

The function M. is the normalized Maxwellian: 

M{y) = -^e^H 2 , (2.3) 
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and the cross section a is symmetric and bounded from below and above by two positive 
constants. The parameter e is intended to go to zero and represents the scaled mean free 

path. The scaled magnetic field is given by f*%, while E = I J is the scaled electric 

field. The gyroperiod is Erf and we shall consider the situation where rj and e go to zero 
simultaneously and the situation where e tends to zero for 77 given then 77 tends to zero. 
Since the magnetic field strongly confines the motion of the particles in the perpendicular 
direction (x,y), we have rescaled this direction with the prefactor rj. This is why the 
orthogonal operator 71 comes with a different scaling from the parallel one T z . By doing 
so, we shall obtain a non trivial motion in the orthogonal direction (guiding center). 

2.1 Scaling 

In this subsection, we shall explain the hypotheses that lead to the scaled Boltzmann 
equation (12. ip . In particular, we highlight the differences of the position scales in the 



parallel and perpendicular directions (with respect to the magnetic field). The starting 
point is the diffusion scaling in the absence of a magnetic field. The unsealed Boltzmann 
equation in this case reads 

d t F + V ■ V X F + ±— ^- • V V F = Q(F) 
m 

where (t, X, V) are the physical time position and velocity variables, while $ is the elec- 
trostatic potential, e is the elementary charge and m the mass of the electron (to fix the 
ideas, we consider electrons). The collision operator Q has the form 

Q(F) = [ £(V, V')[M T (V)F(V') - M T (V')F(V)] dV, 
X is the transition rate and M? is the Maxwellian at the temperature T 
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M T (V) = ( — )- d/2 exp 



m 2K B T 

Kb being the Boltzmann constant. The thermal energy and velocity are respectively 
defined by 

U th = K B T, V th = \ —■ 
V m 

The electrostatic energy is assumed to be of the order of the thermal energy and varies 
over a macroscopic length scale L : 

q$(X) = U th <f>(j). 

The transition rate E is the inverse of a time. The typical value of this time, the scattering 
time, is denoted by t s . So we write 

1 V V 

nvy)==-^—,—) 

r s Vth Vth 

where a is its scaled version assumed to be of order 1. The mean free path L s will be 
defined as the distance that a thermal electron crosses during the scattering time r s 

L s = V th T s . 

This distance is assumed to be much smaller than the macroscopic distance L and we 

shall denote it by 

L s 

The only thing left to know is the macroscopic time scale. Since the kernel of the collision 
operator is generated by the Maxwellian which carries no flux, one has to go to the 
diffusion scaling and define the macroscopic time r m by 
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Using (r m , L, V t h, U t h) as units for time, position, velocity and energy, the Boltzman equa- 
tion can be written 

d t f + -(v ■ V P / - V r • V v /) = 5<£, 

where Q is the scaled Boltzmann operator (12. 2p . 

Let us now look at the effect of a large magnetic field B = Be z . The scaled Boltzmann 
equation contains the additional term — ^(V x B) ■ VyF. An electron submitted to this 
magnetic field spins around it following the equation 

dV t/ n 

m— — = —qv x B. 

dt H 

The period of this precession, called the cyclotron period is given by 

m 

We shall assume that this period is much smaller than the scattering time and we write 

T c = T] 2 T S . 

Of course the motion along the magnetic field stays unchanged and we shall not change 
the length scale in this direction. In the orthogonal direction however, the combined 
action of the electric field and the magnetic field results in the drift of the guiding center 
(the electron motion is a rotation around a center which is now drifting under the action 
of the electric field). The velocity of this drift is given by 

Its order of magnitude is then given by 

v c B , 

wher E± is the perpendicular component of the electric field. 

The last hypothesis that we make is the potential $ varies in the orthogonal direction 
on a lengthscale L± which is not equal to the parallel lengthscale L. Moreover, L± is such 
that the guiding center moves during the macroscopic time by a distance of the order of 
L± itself. In other words, we want 

t? Uth T/ T 

qL± 

Replacing the thermal energy by its expression given above in this section, we find after 
straightforward computations that 

r2 _ T mUth _ 2r2 

L± - ~w ~ v 

Now rescaling the orthogonal position variable by L±, we obtain the Boltzmann equation 

(EH). 



Remark 2.1. From the scaling hypotheses, we see that between two successive collisions, 
the electron precesses a large number of times around the magnetic field. Therefore, one 
might expect that the collision cross section will be averaged along these precessions. On 
the other hand, since the orthogonal lengthscale has been designed in order to see the 
guiding center drift, the limiting equation should exhibit this term. Finally, since the 
parallel motion to the magnetic field does not feel the magnetic field, one is entitled to 
expect a non vanishing diffusion in this direction. This is exactly what we shall prove 
in this paper : drift in the orthogonal direction, finite diffusion in the parallel direction 
governed by a collision operator averaged over precessions around the magnetic field. 



2.2 Notations 

The analysis of the limit e — $■ naturally leads to the study of the operator 

Q r ' = Q- ^^ • V v , (2.4) 

while the limit 77 — >■ 0, involves the cylindrical averaging around the vector e z . To this 
aim, we denote by 1Z(t) the rotation around e z with angle r which is the multiplication 
operator by the following matrix 

(cosr sinr \ 
-sinr cosr . (2.5) 

1/ 

The cylindrical average of a function h = h(v) is defined by 

A(h)(v) :=h(v) = ^- f h(TZ(r)v)dr. (2.6) 

2tt J 



It is clear that A is the orthogonal projector (in the L 2 sense) on the set of cylindrically 
invariant functions. We also define the partial average 

A T (h)(v) : = h r (v) = J- f h(K(s)v) ds. (2.7) 

^ Jo 

We shall moreover define the averaged function a and the symmetric average a by 

_ / 1 f 27r 

°"( v > v ') = 77- / a 0^( T ) V > V ') dr 
2tt Jo 

V (v, v') = J^ p J** a (K(r)v, K(t')V) drdr' = ±- J** a (v, K(r f ) V) dr', 
and the corresponding collision operators 

W)(v) = / a(vy)[M(v)f(v')-M(v')f(v)]dV, (2.8) 

W)(v) = / ^(v,vO[M(v)/(v') -M(V)f(v)]dV. (2.9) 

Jr 3 
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We define the functional spaces 

f = f/\/i / / 

M(y) 



L 2 M--\f = /(v) / / l 4vr^dv < +00 \, (2.10) 



and its cylindrically symmetric subspace 

L 2 M = {fEL 2 M / /(W(t)v) = /(v), VrG[0,27r]}. (2.11) 

We also define the space 

U 2 M = {feL 2 Ml such that, d a f e L 2 M , for |a| < 2} (2.12) 

We shall make the following hypotheses 

Assumption 1. The cross-section a belongs to W 2,OQ (M. 6 ) and is supposed to be symmetric 
and bounded from above and below: 

3ai, a 2 > 0, < cti < a(v, v') = <r(v', v) < a 2 , V(v, v') e M 6 . (2.13) 

Assumption 2. T/iere exzsis a potential V(t, r) m C 1 (1R + ; W /1 '°°(1R 3 )), suc/i t/iat E = 

-V r V. 

We finally define the space 

Mr = {/ = /(', v) / £ Jg^y^ < + oo} . (2.14) 

2.3 Main results 

The main results of this paper are summarized in the following two theorems. The first 
one deals with the limit e — > while rj > is kept fixed. 

Theorem 2.1. Let rj > be fixed. Let T G R+ and assume that the initial data of 
(12. ip . /o, belongs to L 2 M . Then, with Assumptions^ and\^ the problem ( 12. 1ft /ias a 
unique weak solution in C°([0,T],Z/L ). Moreover, the sequence (f £V ) £ converges weakly 
to p v (t, r)Al(v) m L oo ((0,T);L^ /[ )) weaA; * where p v satisfies 

d tPv - div(D"(Vp^ - p„E)) = (2.15) 

and t/ie diffusion matrix 3 V is given by the formula 

X 7 ' = ( ) &em# i/ie on/y solution of \3.7\) . 



Remark 2.2. In the relaxation time case <r(v, v') = -, the matrix B> v can be computed 
explicitly: 



T 



T^+?7 4 



T 2 +r/ 4 

-^- 



T 2 +r) 4 T 2 +r) 4 

V 1) 

In the limit 77 — >• 0, the upper 2x2 block of the matrix reduces to the antisymmetric matrix 

/ 1\ 

X = , whereas its symmetric part is of order rf. The following proposition 

V-i Oj 

shows that these features are still valid in the case of a non constant scattering section o . 



Proposition 2.1. Writing X^ 
i] tends to zero. 



X T ' 

v ^ 1 , we have the following expansions in "H^ as 



x^xf + ^xw + o^ 4 ), 
x^xf + ^x^ + ofa 4 ). 

The diffusion matrix V> v is a definite positive matrix. Its symmetric and antisymmetric 
parts, D^ and D^, have respectively the following expansions 



W s 



as 



v ± ®X ( £ ) dv 
D z 

( 1 0\ 
-10 

V 00/ 



rf^l 



V Bfl + 0( V ') 



o(v 4 



where D z = J v z X z dv and 

TTh(0) - 







X 



(0) 



•2 



X 



(i) 



± 






v z 



dv 



dv 



X 

( 



Vz 



(0) 



X 





(1) 







dv 



dv. 



Remark 2.3. The above theorem shows that the diffusion is of order 0(1) in the parallel 
direction while it scales like rf in the orthogonal one, exactly like in the relaxation time 
coefficient. The antisymmetric part of the diffusion matrix acts essentially on the orthog- 
onal direction and leads to the guiding center motion. However, the above formula shows 
a correction of order 77 of the guiding center motion, which is due to collisions. Indeed, it 
is readily seen that for any three dimensional vector Z 



where 



u' 



>2.Z = u" x Z + Ofa 4 

( -r j jX ( y 0) v z dv\ 
rj J Xi°'v z dv 

1 / 



and (Xx ,Xy ) are the components of Xj_ . This leads to 

-V • (DL(V r p, - a,E)) = V • ( Pri E x u") + 0( V 4 ). 



Theorem 2.2. Under the same hypotheses as for Theorem \2.1\ we now let e and rj 
simultaneously and independently tend to zero (there is no assumption on the relative 
scale between e and r\). In this case, the sequence {f Er, ) erj of weak solutions of (12. ip 
converges weakly to p(t, r)A4(v) in L°°([0,T], L 2 Mv ) weak *, where p is the solution of 



dtp + d z J z + V • (pE x e z ) = 
p(0,r)= / / (r,v)dv 

JR3 



(2.16) 



The parallel current density, J z is given by 

J 2 = -D z (d z p- E z p) 
where D z is the diffusion constant given by 

D z = [ XjVdv, (2.17) 

Jr 3 

and Xi is the zero-th order term of X^ (defined by (14.71) ). 

3 Analysis of the operator Q v . 

Let us consider the space L 2 M introduced in (12. 10J) and define the scalar product on this 
space by 

(f,g)M = / TT dv - 



Let D(Q^) be defined by 

D = D(Q r >) := {/ E L 2 M / (v x e z ) ■ V v / G L 2 M }. (3.1) 

We have the following result 

Proposition 3.1. The operator Q v given by (12.41) with domain D defined by Ii3. 1\) satisfies 
the following properties 

1. For any f G D, we have / Q v (f)dv = (mass conservation). 

</R 3 

2. (—Q V ,D) is a maximal monotone operator. 

3. The kernel of Q v is the real line spanned by the Maxwellian : 

Ker(Q v ) = {f G L 2 M /3a G R such that /(v) = aM(v)}. 



4- Let V be the orthogonal projection on KerQ v . The following coercivity inequality 
holds for any function f G D 

~(Q n (f)J)M>a 1 \\f-V(f)\\ 2 M . (3.2) 

where a± > is the lower bound of a (I2.13P . 

5. The range of Q v , denoted by Im(Q v ), is the set of functions g E L 2 M satisfying the 
solvability condition: I g(v)dv = 0. In addition, for all g E Im(Q r ') there exists 

J ft? 

a unique function f E D satisfying both Q r, {f) = g and / f(v)dv = 0. 

JR3 



Proof. First of all, we recall that the operator Q satisfies all items of Proposition 13.11 (by 
substituting D by L 2 M ). Now Items 1., 3. and 4. of the Proposition follow immediately 
from these properties since we have f RS {e z x v ) • V v (/)dv = J R3 div v [(e^ x v)/]dv = 
and ((e 2 x v) ■ V v f, j)m — 0. In order to prove Item 2., we need to prove that / — Q n 
is one-to-one from D to L 2 M . In order to do so, we rewrite the equation / — Q n (f) = g 
under the following form 

■i(v x e z ) • V v / + (1 + v(v))f = Q + (f) + g, 

where we have written 

Q(f) = Q + (f)-vf, (3.3) 



Q + (f) = ([ a(v,v')/(v')rfv')M(v) 

Jm.3 



(3.4) 

i/(v)=/ cr(v,v')M(v')dv'. (3.5) 

The characteristics of the equation are defined by 

dv 1 . . 

Tt =^ vxe ^ 

v|t= = v 

whose solution is v(t) = 1Z(— ^)v . The solution of the equation can then be defined by 
integration over the characteristics and finally leads to / = J z {Q + {f) + g) where 

n 2 f 2lT 
jr(/i)( v ) =-^— / h{n{T)^^ {l+v{n{s)v))ds dr, 

B r,{v) Jo 

B v (v)=e r ' 2 ^(^(ns)v))ds_ h 

A fixed point argument leads to the existence and uniqueness of / E L 2 M solving the 
equation / — Q v (f) = g in the distributional sense. Since Q is a bounded operator on 
L 2 Ml we immediately obtain that (v x e z ) ■ V v / G L 2 M so that the constructed solution 
is in the domain D of the operator Q v . Let us now prove item 5. Consider g E L 2 M such 

10 



that / g(v)dv = 0. From the maximal monotonicity of Q 77 , for each A > 0, Xld — Q v is 

Jm. 3 
surjective. Therefore, for all A > 0, there exists f\ G D(Q V ) satisfying 

Xh-Q r 'Ux) = g. (3.6) 

Integrating this equation with respect to v gives A / f\dv = / gdv = and with (j3.2p 

we get —{Q r, (f\),f\)M > c^iII/aIIm- Taking the scalar product of (13.61) with f\ in L 2 M) 
one obtains 

A 1 1 /a 1 1 A4 + «i||/a||x < IMI.M- 

This implies that (/a)a>o is bounded in L 2 M . Up to an extraction of a subsequence, there 
exists / G L 2 M such that (f\) converges weakly to /. By passing to the limit A — > in 
the weak form of equation f)3.6p . we get 

1 

V 



2K vxe z )-V v f-Q(f)=g in D'(R 3 ) 



In addition, Q(f) and g belong to L 2 M , then (v x e 2 ) ■ V v / G L 2 M and then / G D(Q V ). Of 
course L 3 fdv = since this is true for f\. Now if there is another solution of —Q(-) = g 
with zero velocity average, the difference of this solution with / is in the kernel of Q v and 
has a zero average. It is necessarily equal to zero. 

□ 

The following corollary is a direct consequence of Proposition 13.11 

/Xi\ 
Corollary 3.1. For all 7] > 0, there exists a unique X r , = ] in (D(Q V )) 3 satisfying 

\ X2 J 

-Q r >(X r ')=\ ^ V± M(v) and [ X"(v)dv = (3.7) 

\ V Z J ^3 

and from the coercivity inequality (13.21) . {X^) v and (r] 2 X. v ± ) v are bounded sequences in L 2 M . 
Namely, we have the following a priori estimates 

WXIWm+W^IWmK—WvMWm. (3.8) 

The following proposition allows to reformulate the problem Q v (f) = —g by using the 
characteristics. 

Proposition 3.2. For all g G Im(Q v ), we have the following equivalence 

- Q v (f) =9**f- L v (Q + (f)) = L v (g) (3.9) 

where Q + is given by (13. 4p and L v is the operator on L 2 M given by the expression: 

L v (f) = -2— [ ^ fCRiT^e 7 ? % v( - n{s)v)da dT, C v (v) = e v 2 C »ins)v)ds _ L ( 3 1Q ) 
^v\ v ) Jo 

11 



Proof. Let (S V ,D(S V )) be the unbounded operator on L 2 M defined by 

D(S») = {/6^/(vxe z ).V v /6^} (3.11) 

S v (f) = ^(vxe 2 ).V v / + Kv)/ V/GD(n (3.12) 

With the decomposition (13.31) . we have Q v = Q + — S v . Proceeding as in the proof of 
Proposition 13.11 (replacing v + 1 by z/), we show that S v is invertible and its inverse is 
nothing but L n . The equation -Q v (f) = g is equivalent to f-(S^y l (Q+(f)) = (S r >)- l (g) 
which concludes the proof. □ 

Lemma 3.1. For all f G Ll oc (M. 3 ), we have: 

\L v (f)\ < ^S - l jf \m(r)v)\e a ^dr. (3.13) 

Proof. This estimate follows immediately from the definition of L r] (13.101) using Assump- 
tion lh D 

Proposition 3.3. Under Assumption [ZJ the solution X* 7 of (13. 7p belongs to 
('H 2 M (M?)) i and X^ and all its derivatives with respect to v of order less than or equal 
2 (up to dividing by the Maxwellian) are polynomially increasing when \v\ goes to +oo. 
Namely, we have 

|X"(v)|<Q?(|v|).M(v) (3-14) 

|^X"(v)|<Q3(|v|).M(v) (3.15) 

3 

for all a = (ai,a 2 ) G N 2 with cti + a 2 < 2 and where d a := 7 d^d" 2 . Here, Q\ and 
Q3 are two polynomials of degrees 1 and 3 respectively (depending on rj). 
Proof. From (13. 9p . we deduce the identity 

X r ' = L v {v r 'M)+L v {Q + {X r ')) > 
where 



1 



v 



'/ - ! V 2 



^7 V ± 



Moreover, we have \Q + (X V )\ < a 2 ( / X v (V)dV)M(v) < c M(v) and with fl3~T3|) one 

obtains (I3.14p . Estimate ( I3.15J) follows by derivation of the above equation with respect 
to v and by using the fact that the cross-section a{y , v') belongs to W 2 '°°(M. 6 ). 

D 

4 Expansion of X 7 ? with respect to 77 

A rigourous expansion of X^ around rj = will be carried out in this section. Corollary 13. II 
provides us with bounds of order 0(1) for X^ and of order 0{\) for X^_. We shall prove 

12 



in this section that X^ is actually bounded in % 2 M . This is due to the fast precession 
around the axis e z generated by the differential term -r(v x e 2 ) ■ V v - Filtering out 
these oscillations is done through the reformulation (13. 9p of the equation satisfied by X. v . 
Moreover, we shall expand X 77 in powers of rf up to the first order. 

Definition 4.1. Let Q be the operator on L 2 M defined by li2.S\) and let Q the one defined 
on L M by Ii2.9\) . Then Q restricted to L M coincides with Q. 

We now list the properties of Q which are completely inherited from those of Q. 



-r2 



Proposition 4.1 (Properties of Q ). The operator Q is a bounded operator on L M 
equipped with the L 2 M scalar product. It satisfies the following properties whose proof are 
immediate and are left to the reader. 



1. For all f G L M , we have Q(f) = Q(f) = Q(f). 

2. —Q is a bounded, symmetric, non-negative operator on L 2 M . 

3. The null set of Q is given by 

JV(<5) = {n.M(v),neM} 

4- Let V : L M — > N(Q) be the orthogonal projection on N(Q), the following coercivity 
inequality holds, 

-m)J)M>a 1 \\f-V(f)\\ 2 M . 

i.e —Q is coercive on NlQ)- 1 . 
5. The range of Q in L M is given by 

Im(3) = NiQ) 1 - = !feL 2 M /jjdv = o\= Im{Q) n L 2 M 

In the sequel, we shall also need the properties of Q considered on the whole space L 2 M 
(whereas Q is restricted to cylindrically symmetric functions). Some of the properties of 
Q are inherited from those of Q but some of these properties like the symmetry are lost. 
The following proposition summarizes these properties. 

Proposition 4.2 (Properties of Q ). The operator Q is a bounded operator. It satisfies 
the following properties 

1. The null set of Q is 

N(Q) = N(Q) = {nM(v),ne R}. 

2. The range of Q in L 2 M is the following set 

Im(Q) = Le L 2 M / J^ ^-dv = o| . (4.1) 
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3. For each g G Im(Q), there exists a unique function f G L M satisfying j fdv = 

lution can be written f = f + 
part offJ = A(f), verifies 



such that Q(f) = g. The solution can be written f — f H — — where the isotropic 



Q(f) = g + Q + ( 9 -= l ), (4.2) 

and where we have denoted Q = Q —V. In addition, there exists a constant C > 
independent of g (and f) such that 



M 



< C\\g\\ M . (4.3) 



Proof. Let us begin with item 1. It is clear that if / is Maxwellian, it is in the kernel of 
Q. Let now / G N(Q). We have Q (/) —Vf = 0. Since the Maxwellian is cylindrically 
symmetric, it is readily seen that Q (f) is always cylindrically symmetric. Therefore, 
/ = =Q (/), is cylindrically symmetric and thus 

W) = QU) = o. 

In view of Proposition I4.1[ / is proportional to the maxwellian. 

In order to prove items 2. and 3., we consider g G Im(Q), and / G L 2 M such that 
Q(f) = g. Decomposing / and g as / = / + hf (resp. g = g + h g ) where hf = h g = 0, we 
have 

Q(7) + Q + (hf)-9 = h 9 +Vh f . 

The left hand side of the above equality is cylindrically symmetric, whereas cylindrical 
averages of the right hand side vanish. Therefore, both are equal to zero. We deduce that 



h, = - 



h g _ 9~ 9 



V V 

whereas / satisfies 

Q(f) = 9-Q^(hf) 

which can be rewritten Q(f) =~g — Q (^=^) and also 

^(7+|) = g + (|). 

This equation has a solution if and only if the right hand side has zero velocity average 
(see Proposition 14.1)) . A straightforward computation shows that this is equivalent to 

the condition f ^ = 0. Let U8 finally 8 how e 8 t„„ate » Remark first that for 

JR3 V _ _ 

any / G L 2 M with f = f + h f and hf = 0, we have \\f\\ 2 M = \\f\\ 2 M + \\h f \\ 2 M and then 
H/H.M ^ II/II.M- Besides, equation (14. 2p is posed on L M then, from the coercivity of — Q 
on N^Q) 1 - (see Proposition 14. ip . there exists a constant C\ > such that 

m < Ci(||3+g^(^)|U0 < Ci(\\g\\ M + \\Q + {^)\\m) < C[\\g\\ M . 
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In addition, since ht = — — — and with (12.1 3D . there exists Ci > such that II /if II m < 

y 

C2\\g\\M an d estimate (I4.3p holds. 

D 

Definition 4.2. We define the gyration operator Q by G(f) = (v x e z ) • V v / ; with domain 
D. 

The kernel of this operator is nothing but the set of cylindrically symmetric functions. 
It is clear that its range is contained in the set of L 2 M functions which have zero cylindrical 
averages. The following lemma shows that both sets coincide. 

Lemma 4.1. Let g G L 2 M with zero cylindrical average (A(g) = 0). Then, the set of 
solutions of the equation 

Q(f)=9, (4-4) 

in nonempty. Any solution f can be written 

f = fi+A l (g) 

where /i si an arbitrary cylindrically symmetric function (which means that A(f\) = f\) 
and the average operator A\ is defined by 

A l (g) = —- g{Jl{r)v)— T dr = — g{K{r)v) v{K{s)v) ds dr. (4.5) 

v\y) Jo 27n/ j j 

Moreover, we have 

A(vAi{g)) = 0. 

Proof. Let us consider the problem 

G(f v )+V 2 vfv = 9- 
This problem has a unique solution f v defined by 

fv = - 2 L v (g) = p^-r J** g{K{T)vy 2 %^ s ^ ds dT 
V L/ v\ v ) Jo 

1 /" 27r 

g(n(T)v)[e n2 & u{n{s)v)ds - l}dr 



C v (v) jo 
where 

C V { V ) = e V 2 .t»(K(s)v)ds _ x = e 2^V _ L 

Passing to the limit rj — > 0, we find that f v converges almost everywhere towards Ai(g). 
Moreover Aiyfvf^) = A(g) = 0, which proves after passing to the limit rj — > that 
A(vAi(g)) = 0. Of course, any other solution of the equation Q(f) = g is obtained by 
adding an arbitrary cylindrically symmetric function. 

□ 

Let us now consider the problem 

QU) = 9 



A(Q(f)) = h. (46) 

It is clear that the above system does not have a solution unless A(g) = and h G L M 
with L 3 hdv = 0. The proposition below shows that these conditions are sufficient and 
gives the general solution of the problem 
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Proposition 4.3. Let g and h be in L 2 M such that A{g) = 0, A(h) = h and J hdw = 0. 
Then the problem \-j.b]j has a unique solution up to a Maxwellian (the Maxwellian is the 
only solution of the homogeneous problem). The unique solution f which has zero average 
is given by 

f = f l+ A l (g) 

QU) = -vA l (g) J rh 

with /i is an arbitrary cylindrically symmetric function. Moreover, the mapping (g,h) h-> 
/ is linear and continuous both on L 2 M and H 2 M , ie. 

H/lk < c{\\h\\* M + NI*J' ll/llflj, ^ ^ll^ll^ + IMI**,)- 

Proof. First of all, Lemma [4.11 shows that / = /i + Al(<7). Besides, we have 

Since /i is cylindrically symmetric, then «4.(f/i) — F/i = 0, which leads, thanks to the 
identity A{vA\{g)) = 0, to Q(f) = —VAi(g) + h. It is now enough, in view of the 
solvability condition (14. Xp to check that j =(—VAi(g) + /i)dv = 0, which is readily seen. 

□ 

4.1 Expansion of X^ 

The idea is to make a Hilbert expansion of X^ in powers of rj 2 : 

i 

We have the following equations 

S(xi 0) ) = 0, 

^(Xi 1 )) = Q(X z 0) )+v z M, 

gix z i+l) ) = g(xi 4) ), z>i, 

where we have taken advantage of the fact that A{v z M) = v z K\. The solvability condi- 
tions become 

A(Q(X^)+v z M) = 

A(Q(X$>)) = 0, i>l. 

Therefore, the terms of the expansion can be successively computed by solving the prob- 
lems 

G(X Z 0) ) = 

^(g(xl 0) )) = - Vz m. 

G{X { P) = Q(X^)+v z M. 

A(Q(X z iy )) = 0. 

g(xi i+1) ) = Q(xi l) ) 

A(Q(X z l+1) )) = 0. 
16 



From Proposition 14. 3[ we deduce that Xi are uniquely determined by the above equations 
and are in TL 2 M . 

Proposition 4.4. Let X* 7 be defined by ([ff.Tj ). The function X\ is uniquely determined 
by 

-^(X^) = v z M, f xW(v)dv = (4.7) 

</R3 

and X z is the unique solution in L 2 M of 
- Q(X«) = uA 1 (Q(X^)+v z M) = ^(Q(Xf ) -5(X(°))), / X«(v)rfv = (4.8) 



where A\ is the averaging operator defined by (14.51) . Then the following expansion holds 

,2 



in the strong topology ofH 2 , 



X« = xW+ v 2 XW+0(ri*). (4.9) 

Proof. We write 

X^Xf+^XW+r^V. 

From the definitions of X z , it is readily seen that 

£(r") = r / 2 Q(r")+Q(X«) 

A(Q(r^)) = 0. 

Hence, we have that 

Q^) = -VA l {r, 2 Q{^) + Q{X^)). 

Using estimate (14. 3p . we obtain for some constant C > independent of r\ 

PII^<^ 2 p|| Lli + ||xW||^) 

which shows the boundedness of r v in L 2 M . Differentiating the above identity with respect 
to the velocity variable, on can show simply the boundedness of r v in V? M . □ 

4.2 Expansion of X^ 

In this subsection, we expand the orthogonal part of X 71 . We follow the same steps as for 
X^ and write the equations satisfied by the formal expansion of X^_. Namely, we have 

e(xf) = v ± m, 

0(x£>) = Q(xf ), 
£(X? +1) ) = Q(X«), i>\. 



17 



As for X^, the above system can be reformulated as follows 

^(Xf) = v ± M 

A(Q(xf)) = 0. 

g^) = Q(xf) 

A(Q(X^)) = 0. 

S(xf :) ) = g(x«) 
^(g(x? +1) )) = o. 

Proposition 4.5. T/ie expansion 

Xl = X^ + r ] 2 X { l ) + 0(r ] i ), 
holds true in rl 2 M . Morever, X\_ satisfies 

Xf - A(xf) = -X(v ± M), 



with 1 



1 
-1 



The proof of this proposition is identical to that of Proposition 14.41 and is skipped. 
The only thing which has to be checked is the expression of the anisotropic part of Xj_ 
which can be deduced from the identities 

-v ± M = g(l(vj_M)), A(I(v ± M)) = 0. 

5 Proof of the main theorems 

We begin this section by making precise the definition of weak solutions of the Boltzmann 
equation and give quite standard a priori estimates on this solution. For the moment, 
e and rj are arbitrary and all the constants are, unless specified, independent of these 
parameters. 

Definition 5.1. Let T e M.* + , a function f £r > e Lj oc ([0,T] x R 6 ) is called weak solution of 
(12.11) if it satisfies 



T 



f 



= '/ 



d t ip H h 



srj 

T 



for all $ e Cl([0,T) x 



dtdrdv 



£ 2 V 2 Jo 



f £v g(i))dtdrdv 



Q(f £r] )ipdtdrdv + / f (x,v)ip(0,x,v)dvdv. [5.1] 
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Theorem 5.1. Assume that fo G L 2 M (drdv). Then, Ve, 77 > andVT G R+, t/iere exists 
a unique weak solution f £r} G (7°([0, T],!/^ (drdv)) 0/ (12.11) . In addition, the following 
mass conservation equation holds 

d t p £V + div r J E " = (5.2) 

where 

p £ \t,x) = f rdv, J? = - f vifdv, J? = — I v±f £V dv, (5.3) 

Jm.3 e j R a srj j R 3 

and there exists a constant C > independent of e and 77 such that 

\\f £V (t)\\ Lliidrdv) <C VtG[0,T], (5.4) 

\\r - p^M\\ LH[onL%t(drdv)) <Ce (5.5) 

An immediate consequence of the integrability properties of the function f £V , with 
respect to v is that the weak formulation (15. ip is satisfied by test functions which are not 
necessarily compactly supported in velocity. More precisely, we have 



Corollary 5.1. Let f er] be the solution exhibited in Theorem \5.1[ Then the weak formu- 
lation 115. 1\) is satisfied by test functions lying in the set 

T — {ip(t, r, v) G C 1 ([0,T) x M 6 ) compactly supported w.r.t. t and r and 

3neN,CeR + , |V t , r , v ^(t,r,v)|<C(l + |v|) n }. (5.6) 



We shall not prove this corollary which is straightforward. The proof of Theorem 15.11 
is also classical (see for example [32]). We show here how the a priori estimates (15. 4p 
and (15.51) can be obtained. We simply multiply the Boltzmann equation (12. ip by j^- and 
integrate with respect to r and v. Straightforward computations show that this leads to 
the entropy inequality 

d f^ drdv- I d t V^-drd\ 



dtj 2M V J 2M V 

1 f f £V rci 

=- Q^nj—dvdv <-^\\r - p £ ^m\\ 2 l2 . 

e z J Jviy e z M v 

Thanks to the boundedness of the potential V and its derivatives, estimate 05. 4p follows 
by applying Gronwall Lemma and then the bound ( I5.5I) follows immediately. 



5.1 Proof of Theorem 12.1 



In this section, we are interested in the limit e — > 0, while 77 is fixed. To this aim, we follow 
the usual approach based on the moment method. Namely, thanks to estimates ( 15. 4p and 
(15. 5p . the sequences f £V , (p £V ), (J^) and (J^ 7 ) are bounded with respect to e. One can 
find p v (t, x) G L 2 ((0, T) x M 3 ) such that p £V converges weakly to p v in L 2 ((0, T) x M 3 ) and 
f £r] ->> p v M in L°°(0, T; L 2 M (R 6 )). Similarly, J e?? converges weakly to ^ in L 2 ((0, T) x : 
We have of course, the mass conservation equation 

d t Pn + div r 3 V = 0. 
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Now the only thing left to show is to identify 3 V . The way to proceed for identifying J v 
is standard (see for instance [U EJ ESI El E2J QUI El El II])- The Boltzmann equation can be 
rewritten 

at r\ e 

Let now X 77 * be the solution of 

qv*(xv*) = - ( ^ v± I M(v) and / X"*(v)dv = 0, 
\ v z ) J& 

where Q v * := (Q v )* is the adjoint in L 2 M of Q v . It is readily seen that Q v * has the same 
expression as Q v except that rf is replaced by — rf. We then have the following result. 

Corollary 5.2. The following expansion holds in l-i 2 M 

X 7 !* 

We multiply (15. 7p by —^— and integrate w.r.t. v. The right hand side of the obtained 
equation is nothing but 

i fQ ri (f eri )xr dy= i [ Q*\xr)r ^_ jev 



M e ./ M 

Therefore 

Thanks to the bounds (15. 4p and (13.141) . we can pass to the limit in the above equation as 
e tends to zero and obtain 

j^T z ( Pri M)dv + -J -^T ± ( Pv M)dv, 
which leads to the expression 

-J? = [j X?v z dv](d zPn + Pv d z V) + -[fx?v ± dv] ■ (V r±Pv + p v V r± V). 

Multiplying the Boltzmann equation (15.71) by ^ and proceeding analogously, we find 
out that 

-Jl = 4 J XYv t dv]{d tPv + Pv d z V) + [Jxf <g> v ± dv](V r± p„ + P v V r± V). 

Using the property f X^*v z dv = f v z X^dw as well as f X^*v±dy = rf f t> z X^_<iv and 
analogous formulae for the terms appearing in the expression of 3\, we finally obtain the 
expression 

JT' = -0 , >(V r p„ + p„V t .V) 
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where 



One can then expand V> v as 



where 



D^ = B v ± + W Z + r}W z± 



5.2 Proof of Proposition 12.11 

Let us expand first W]_. Following the expansion of X^_* and that of X^_, we obtain that 



where 



and 



»r=-ff x fUf;)^/(M 8 (f ,* 



»«. //x<'MJv,\ dv= f(M /x« |(Jv _ 



o / \ o ) ~ \ o 



It is clear that H> ± is antisymmetric while Dj_ is symmetric and we have 

.<»> = _ K Xf -4(X«°') ^ f vx , dv . 




Since X^ - ^(xf ) = -Z(v ± .M), we have 



D ( , 0) 



/ l(v ± M)®v±dv | ( X- I (v ± ®v ± )Mdv \ 
° / \ ° ° / 



1 

-10 



where we have noticed that / vj_ <g> v^Aidv is nothing but the identity matrix on 
Let us now expand D^. We have 

W z = Df + V 2 3^ + Q{rf) 
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■•/l 



where 

/ \ 

b>(°) = D z 

\o 1/ 

which D 2 = J Vz'X.i dv. Comparing the expansions of X^and X%*, we deduce that W z 
is equal to zero. Let us now expand D^ ± . Of course we have 

D?x = »S + -M'i + 0('f 4 ), 



#._ /fowxrw/fxrw 



and 



where 

i o) \ dv _ r( xf 

We can now make explicit the expansion of D^ = D^ + D^ s where the indices s and as 
stand for the symmetric and antisymmetric parts 



ri 2 fv ± ®X i l ) dv 

D ? 



ZmW _l /nr„4\ 



+ n a w;l + o(ri 



5.3 Proof of Theorem [2721 

Lemma 5.1. The density p £V and the parallel current J z v (15.31) are bounded in L 2 X ([0, T] x 
R 3 ) with respect to e and r\. In addition, there exist p and J z in L 2 ([0,T] x IR 3 ) such that, 
up to extraction of subsequences, we have 

p eT] — ^ p weakly in L 2 X 

fi^pM inL°°([0,T],L 2 M (drdv)) 
Jf-J, inL\ x . 

However, estimates (15. 5p and ( I5.4J1 do not give a bound to J^ and are not sufficient 
to close the limit equation or to find relations between the corresponding limits of the 
current and the density. To deal with this lack of compactness, we have to filter out the 
oscillations generated by the magnetic field in the orthogonal direction. 

Proposition 5.1. Let J^ be the perpendicular part of the current given by fl5.3|) . Then, 

J^-^Jl mD'([0,T)xl 3 ) 

with 

J± = (pE ± - V r± p) x e z 
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rpC* 
Proof. Multiplying (12. ip by ^ and integrating with respect to v, we get 



-J 



:>t 



7 R 3 A* ,/ R 3 M '7r3 -M 



nX^* 
More precisely, by taking a test function — —^(fi(t,x) in the weak formulation (15. ip . 

where </>(£, x) G C£([0,T) x R 3 ) we have 

T /" 

/ J±(t,x)(J)(t,x)dtdx 

o </R 3 

= £7 1 f £V -j-d t <f)dtdrdv+ / / f er >T±_{—tT-)dtdrdv 

Jo Jr6 -M Jo Jr6 M 

+V f f f £r, T z (^)dtdrdv + e V [ f Q {x,v)^±-<l>{0,x)dxdv. 

JO JR6 M J R 6 M 

In view of Corollary 15.21 X*_[* — > — ~X. ± strongly in H 2 M while (f £V ) converges to pAi 
weakly L 2 (0, T, L 2 M {drd\)). Therefore, we can pass to the limit in the above formula and 
get 

r T r r T r X^</> 

lim / / J £ ^(t,x)(p(t,x)dtdx = - / pMT±(^ L —)drdvdt 

erf-* Jo Jr3 J J R e M 

= pJ(V r± 0-V r± y0)drdt, 

Jo Jr 3 

where 

X=[ ) =- [ xfsvidv. 

\ -1 / J^ 3 

This concludes the proof of the proposition. 



D 



The only thing left to do is to give the expression of the limiting parallel current. To 
this aim, we use the test function —rr<f>(t, x) in the weak formulation where is a regular 
compactly supported function 

/ J £ z v (t,x)0(t,x)dtdx 

ii Jr 3 

/ rT z {-^r)dtdrdv + e f (x,v)-f-<f ) (0,x)dxdY. 
o Jr & JV[ Jr 6 JV[ 

rT p ~YV* 1 pT f X 7 ?*^ 

+£ / / f £n ^rd t (t)dtdTdw+- / f £7 iT±(^-r)dtdrdv. 
Jo Jrb M v Jo Jr6 M 

(5.8) 
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The first term of the right hand side converges as e and 77 tend to zero towards 
/ pMT z (^-y-)dtdrdv = / / pD z {d z (j) - (j)d z V) dtdr. 

Jr 6 M JO Jr 3 

The second and third terms obviously tend to zero. The following lemma shows that the 
last term of ( 15. 8 p converges to zero. 



Lemma 5.2. Let <f>(t, r) be a compactly supported C 2 function. Then, we have 

e,»J->0 T] J J R6 



lim - / / f £r <TA^-r)dtdvdv = 
^o?l Jo Jr« M 



Proof. Since we have XIT = Xi + 0(rf) in H^, strong, we can replace X!f by X« . Now, 
we remark that Xi is cylindrically symmetric so that 

Therefore, one can find a test function -0(t, x, v) in the set T such that 

The fact that ^ lies in T comes from the formula ip = Ai(T±( jb )). Using this function 
-0 in the weak formulation, we obtain 

i r r /• _ ,X (0) 



^ Jo Jm 6 



x ch \ r r 

f r 'TA- I rT-)dtdYdw = - / f £1 iQN,)dtdrdv. 

v -M ^Jo </R6 



Using the weak formulation (15.11) . the right hand side of the above identity can be imme- 
diately estimated as 0(e), which tends to zero as e and rj tend to zero. 

□ 

6 Concluding remarks 

We have proven in this paper that the diffusion limit of the Boltzmann equation with 
an ultrastrong magnetic field leads to a diffusion equation in the parallel direction and 
a guiding center motion in the orthogonal one. The proof has been done for the linear 
Boltzmann equation and in the case of a constant magnetic field and involved the analysis 
of the operator involving collisions and gyrations around the electrostatic field. If the 
magnetic field has a constant direction but smoothly varies in position and time while 
staying away from zero, the analysis can be carried out without difficulty and the results 
can be generalized. A more difficult problem appear if one couples the Boltzmann equation 
to the scaled Poisson equation or when the collision operator has nonlinear features. In 
this case, the method used in this paper might not be sufficient and the use of double 
scale limits might be necessary as it is for gyrokinetic limits. 
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